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INTRODUCTION 

The behavior of QCD at finite temperature is of great 
theoretical and experimental interest. Traditionally, fi- 
nite temperature calculations of a field theory are per- 
formed in Euclidean space with periodic boundary con- 
ditions imposed in the time direction. However, the same 
information can be extracted from a Minkowski space cal- 
culation by examining the density of states in the spec- 
trum as a function of energy. In particular, an exponen- 
tially increasing density of states implies a finite critical 
temperature Tc . There is some experimental evidence 
of this exponential increase in the hadron spectrum 2]. 
Recently, there has been renewed interest in the lightcone 
formulation of thermal field theory in Minkowski space 
I^J. In this letter we investigate the aspects of a light- 
cone formulation of QCD relevant for finite temperature 
and, in particular, its implications for the string picture 
of QCD. 

To leading order of the l/Nc expansion of QCD 0, 
we compute the density of glucball states per unit mass 
via the transverse lattice approach to non-abelian gauge 
theory This is a Hamiltonian method that combines 
light-front quantization in the two "longitudinal" space- 
time directions, = ±x^)/V^, with a lattice in the 
remaining "transverse" directions, {x^,x^}. It has been 
successfully used to compute the properties of the light- 
est glueballs and mesons (for a review, see Ref. 0). In 
these calculations, effective Hamiltonians on coarse lat- 
tices were tuned to optimize Lorentz covariance of the 
lightest boundstates. The method also yields the spec- 
trum of heavier boundstates. A rapid exponential rise in 
the density of states suggests that there may be enough 
states at relatively low masses, where the coarse lattice 
calculation is still accurate for aggregate quantities, to 
observe this behaviour unambiguously. 

Using the detailed methods described in Refs. 0, we 
demonstrate that in an intermediate mass range, ~ 1 — 3 
times the lightest glueball mass, the density of states 
does indeed show a degree of universality with respect 



to cutoffs. We find an exponentially increasing density 
of states, implying a finite-temperature Hagedorn transi- 
tion. Although the method is not well-suited to produc- 
ing an accurate prediction of the critical temperature, we 
make a numerical estimate and produce analytic bounds 
on it. One application of interest is to the string pic- 
ture of large- A'^c gauge theory. Results from previous Eu- 
clidean lattice computations indicate a Tc lower than 
that of the Nambu-Goto string model 9j, suggesting that 
a QCD string in the large- A^c limit has more worldsheet 
degrees of freedom than this free bosonic string. Our 
calculations suggest that this is due at least in part to 
longitudinal oscillations. This is also consistent with the 
Nambu-Goto string description appearing in a distinct, 
unstable phase of transverse lattice gauge theory discov- 
ered by Klebanov and Susskind [lfij |. 



TRANSVERSE LATTICE GAUGE THEORY 

The transverse lattice is formulated in the following 
manner. Two coordinates x", a € {0,3}, are continu- 
ous while two directions x = {x^,x^} are discretized as 
a lattice, spacing a. The longitudinal continuum gauge 
potentials A°'{x.) lie at sites x. The transverse fiux link 
fields Mr(x), r g {1,2}, lie on the link between x and 
X -I- ar, representing the gauge fields polarized in the x 
directions. In general there are also fermi fields ^'(x), but 
we will not need these for the glueball spectrum at large- 
Nc- For a coarse transverse lattice, the strategy is to 
perform a color-dieletric expansion [llj of the most gen- 
eral lightcone gauge Hamiltonian, renormalisable with re- 
spect to the continuum coordinates a;", in powers of M^-. 
Provided these link fields are chosen sufficiently heavy, 
one can truncate this expansion to study the low lying 
hadron boundstates dominated by just a few particles 
of these fields. The remaining couplings in the effective 
Hamiltonian can then be accurately constrained by op- 
timizing symmetries broken by the cutoffs in low energy 
observables. Typically this method is viable in a window 
between small lattice spacings, where the fields become 
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too light to justify the color-dieletric expansion, and large 
lattice spacings where Lorentz covariance breaks down 
uncontrollably. 

For completeness, we give the Lagrangian, used in 
Ref. 0, that we shall employ here. It contains all al- 
lowed terms up to order (M)^ for the large- A^c limit: 
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where F"'^(x) is the continuum field strength in the 
{x°,x^) planes at each site x. We have defined Air = 
M^r and hold G^Nc finite as Nc oo. In the iV^ ^ oo 
limit, glueball states consists of connected loops of trans- 
verse flux links and are absolutely stable to decay. 

The reader is referred to Ref. j2] for details of the con- 
struction of the lightcone Hamiltonian and Fock space, 
the renormalisation, and the determination of various 
couplings appearing in Eqn. These couplings are 

accurately constrained by optimizing covariance of low- 
lying glueball eigenfunctions and rotational invariance of 
the heavy-source potential. The relatively large num- 
ber of couplings means the Hamiltonian is highly 'im- 
proved' and can give cut-off independent results on quite 
coarse lattices. Taking the fundamental scale to be the 
string tension a, Refs. j3| investigated individual glueball 
masses up to about 1.5 times the lightest glueball mass 
and the method was shown to produce accurate results, 
even for lattice spacings of order I/Vct. We investigate 
here whether any part of the higher excited spectrum, 
probed at finite temperature, is similarly under control. 



DENSITY OF STATES 

We employ a DLCQ cutoff K in the direction 
to obtain a finite Fock space basis at finite transverse lat- 
tice spacing a. Therefore, we must study the stability of 
results under variation of both K and a. Figs, d and El 
display the distribution of glueball masses Aft, in units of 
the string tension for the t**^ state in the spectrum. 



The results in Fig are shown for fixed a = 1.44/v^, 
the lattice spacing at which the best covariance of the 
low- lying spectrum was obtained for Lagrangian Q , and 
increasing values of the DLCQ cutoff K ^ 8, 10, 12, where 
_ftr — > oo is the limit of infinite dimensional Fock space. 
The largest basis at iiT = 12 contains 654,948 states in 
total, though only the lowest 20,000 are plotted for file 
storage reasons. At the largest masses displayed, the 
DLCQ cutoff is introducing non-universal artifacts be- 
cause it cuts off the maximum number of Af-particles in 
a state. At lower masses however, there is universal be- 
havior of the slope of the distribution. We repeated the 
calculation for K = 12 aX four different lattice spacings, 
shown in FigEl the couplings of ^ re-optimized for co- 
variance of the low-lying spectrum at each. Again, for the 
range of masses plotted, the slope is relatively universal. 

The approximately universal linear growth observed 
in this region implies that the density of states per unit 
mass, p{M) ~ dt/dM, grows exponentially with mass M 



p{M) 



(2) 



This exponential behaviour establishes itself already in 
the first 100 or so states, occuring at less than twice the 
lightest glueball mass, which is within the expected range 
of validity of the method 01 . It is perhaps surprising from 
the figures that the universal growth appears to extend 
to much higher masses. One would expect the individual 
highly excited glueball masses to shift around relative to 
their true values due to cutoff dependence. However, an 
aggregate quantity like the density of states appears to 
be less sensitive to cutoff dependence. 

CRITICAL TEMPERATURE AND STRINGS 

It is well known 1] that exponential growth of the 
density of states leads to a diverging canonical partition 
function for temperature T > Tc- To determine an ac- 
curate value for Tc, one needs to know the density for 
asymptotically high masses M, when power corrections 
to the expression ((SJ at lower masses can be neglected. In 
principle, knowledge of the functional form of power cor- 
rections could enable one to extract accurately from 
lower masses, but in practice the form is unknown for 
QCD. Therefore, although exponential growth is clearly 
occuring, we can only extract a rough estimate for 
from the numerical data, based on linear fits to FigsQ] 
and El 

There are three systematic sources of error in our esti- 
mate. The first two come from the variation of the data 
in Figs, m and El which we take as indicative of the error 
due to residual dependence on the cutoffs K and a. The 
third comes from the linear fits themselves, since power 
corrections to Eqn. (0) in the actual data will make the fit 
result for Tc depend upon the range of masses chosen for 
the fit. Those intervals are naturally limited above and 
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FIG. 1; ln(t) versus the mass Mt of the t"^ glueball in units of the string tension, for K — 8, 10, 12 and a^/a — 1.44. Lighter 
data correspond to larger K. 
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FIG. 2: ln(t) versus the mass Mt of the t**^ glueball in units of the string tension, for a-Ja — 1.09, 1.21, 1.33, 1.44, and K — 12. 
Lighter data correspond to larger a. 



below: for all the data, smooth linear behaviour does not 
start until the block of points beginning at M k, Q^/a\ 
DLCQ artifacts limit the maximum useful value of i, 
roughly Int = 6, 7,8 for the if = 8, 10, 12 data respec- 
tively. Within these limits, we performed linear fits to 
the data sets in the figures for a variety of different mass 
intervals. The maximum variation of Tc with interval 
chosen was of order ±15%, which we take as indicative 



of the error made by using simplified fitting function 
As a result, we estimate Tc — 0.63(1)(9)(9)-\/ct, where the 
first error is due to K ^ second due to a, and third from 
the fits. 

We can also deduce analytic bounds on Tc. In Ref. [llj |. 
we derived in transverse lattice gauge theory an approxi- 
mate analytic formula for the number dn of highly excited 
glueball states of mass squared Af^ = 27r(Tn, n integer. 
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The derivation was based on the assumption that the 
glueball masses were governed only by longitudinal dy- 
namics of the transverse link fields and a fixed number of 
such links in a state. Assuming a transverse degeneracy 
factor (jP for the number of distinct states with p links, 
the formula of Ref.^Jl trivially generalises to 

dn = exp{2^-nLi2[-.g2]} , (3) 

where Li2 is the dilogarithm function. This leads to a 
Hagedorn temperature 

The maximum degeneracy one could expect is from a 
long free chain of p links, which on a square transverse 
lattice gives 3 = 4. This provides an upper bound on the 
density of states and hence a lower bound Tc — 0.5A^/a. 
Conversely, setting g = 1 corresponds to neglecting trans- 
verse structure altogether. This is a lower bound on 
the density of states, implying Tc = 1.375-^0^ is an up- 
per bound. The numerical solution investigated above is 
somewhere between these extremes; in reality there are 
non-trivial longitudinal and transverse dynamics. 

Recent large-A'c extrapolations of Euclidean lattice 
computations ;;8| are more suited to providing accu- 
rate determinations of the critical temperature, obtain- 
ing Tc — 0.596(4)y^. This is less than the result, 
Tc/y/o' = v^ZIj'^ ~ 0.69, for the non-interacting Nambu- 
Goto string This would suggest that the large- A^'c 

gauge theory flux string has more worldsheet degrees of 
freedom than a simple bosonic one, which has only trans- 
verse oscillations. In making this comparison based on Tc 
with a free string, it is important to use results from the 
large- TVc limit, since only in this limit do lightconc flux 
strings not split and join. At finite Nc, by examining only 
Tc one cannot separate physically distinct world-sheet ef- 
fects from string interaction effects proportional to 
Also, the Hagedorn determination of Tc directly from the 
density of states is strictly correct only for a gas of free 
hadrons / strings . 

Although Euclidean lattice results provide an accurate 
value for Tc, they do not explain why it differs from the 
bosonic string value. Transverse lattice gauge theory can 
provide an explanation with the help of an older analysis 
by Klebanov and Susskind They noted that choos- 
ing a sufhciently tachyonic link field mass fj,h (see EqnQ 
would drive the theory into a new phase, where glueball 
states were unstable to becoming infinitely long chains 
of transverse links, each link with frozen longitudinal dy- 
namics. In this phase, there are only transverse dynamics 



(this is the complete opposite of the scenario used above 
to derive the bounds) . The spectrum of this phase is ex- 
actly that of the Nambu-Goto string |lCl Il4l|. Thus, at 
least some of the extra degrees of freedom present in the 
QCD phase of transverse lattice gauge theory, which can 
lower the critical temperature, appear to be longitudinal 
oscillations of the flux string. 
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